BCS-BEC crossover in a quasi- two-dimensional Fermi gas 
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We consider a two-component gas of fermionic atoms confined to a quasi-two-dimensional (quasi- 
2D) geometry by a harmonic trapping potential in the transverse direction. We construct a mean 
field theory of the BCS-BEC crossover at zero temperature that allows us to extrapolate to an infinite 
number of transverse harmonic oscillator levels. In the extreme BEC limit, where the confinement 
length exceeds the dimer size, we recover 3D dimers confined to 2D with weak repulsive interactions. 
However, even when the interactions are weak and the Fermi energy is less than the confinement 
frequency, we find that the higher transverse levels can substantially modify fermion pairing. We 
argue that recent experiments on pairing in quasi-2D Fermi gases [Y. Zhang et al., Phys. Rev. Lett. 
108, 235302 (2012)] have already observed the effects of higher transverse levels. 



Recent experimental advances have made it possible 
to study quasi-two-dimensional (quasi-2D) atomic Fermi 
gases in a very controlled manner [1-8J. Such simple 
quasi-2D systems may provide useful insights into struc- 
turally complicated unconventional superconductors like 
the cuprates, where superconductivity originates in the 
copper oxide planes [9;. Quasi-2D Fermi gases are also 
of fundamental interest, since they are the marginal case 
of the Mcrmin- Wagner theorem and thus have modified 
superfluid properties [TQl [Tl] . In addition, the quasi-2D 
geometry can strongly affect fermion pairing within the 
superfluid, as we investigate in this paper. 

In cold-gas experiments, atoms may be confined to one 
or more quasi-2D "pancake" structures using a ID optical 
lattice. The system can then be tuned from 3D to 2D by 
increasing the confining lattice potential. For sufficiently 
strong lattices, the confining potential for a single quasi- 
2D layer can be modelled as a harmonic oscillator poten- 
tial, V(z) = ^mui 2 z 2 , in the transverse z direction, where 
to is the atom mass. When the temperature k B T <C fkv z 
and the Fermi energy sf <C huj z , the atoms will reside 
in the lowest harmonic oscillator level in the absence of 
interactions, and the gas is considered to be kinemat- 
ically 2D. The short-range interatomic interactions are 
also renormalized to yield an effective 2D s-wave scatter- 
ing amplitude and associated two-body bound state with 
a binding energy e b that depends on both the 3D scatter- 
ing length and the confinement [12 H2]. Thus, by varying 
the ratio Eb/ef, one can explore the crossover in 2D from 
weak BCS pairing (eb/^f *C 1) to the Bose Einstein con- 
densation (BEC) of dimers (e b /e f > 1) [HI HI]. How- 
ever, the interactions also mix in higher harmonic levels: 
for instance, in the BEC limit, dimers will be smaller than 
the confinement length l z = y/KJmuj~ z once £b > f^z, so 
that they essentially become 3D bosons confined to quasi- 
2D. Here, we are interested in how the confinement can 
impact pairing and lead to a departure from 2D behavior 
throughout the BCS-BEC crossover. 

We focus on zero temperature, where there is a well- 
defined condensate in 2D, and we construct a mean- field 



theory that generalizes the 2D results of Randeria et 
al. [2] to quasi-2D. We expect the mean-field approxi- 
mation to be reasonable since it appears to be consistent 
with recent experiments in the 2D limit (ef <S hui z ) [7]. 
The BCS regime of the quasi-2D Fermi gas was previ- 
ously studied in Ref. [TB] using a mean-field Bogoliubov- 
de Gennes calculation that included the lowest three har- 
monic oscillator levels. Our approach, however, allows 
us to extrapolate to an infinite number of levels and 
thus explore the entire BCS-BEC crossover. In the limit 
e b 3> huj z , we find that our calculation recovers weakly 
repulsive 3D bosons confined to quasi-2D. However, even 
for weak interactions, eb *C huj z , we find that higher har- 
monic levels can substantially modify fermion pairing as 
we perturb away from pure 2D and Ef approaches Huj z . 
We determine the radio frequency (RF) spectra for the 
quasi-2D Fermi gas and show that recent measurements 
of pairing in the quasi-2D regime Ef ~ hu z [H] are con- 
sistent with effects due to higher transverse levels. 

In the following, we consider a two-component (f, \) 
Fermi gas interacting close to a broad s-wave Feshbach 
resonance in 3D. Under harmonic confinment in the z- 
direction, the motion of each atom can be parameterized 
by its 2D momentum k in the x-y plane and the harmonic 
oscillator quantum number n in the transverse direction. 
The many-body Hamiltonian is thus (setting the system 
volume to 1): 

H = ^ ( £kn ~ M)cL<T C kr«r (1) 
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where ek« = k 2 /2m + nui z are the single particle energies 
relative to the zero-point energy of the n = state (we 
now set H = 1), and /i is the chemical potential. Note 
that we assume the mass and chemical potential (and 
thus the particle density) are the same for each spin a. 
Since the short-range interactions only depend on the 
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relative motion, we obtain the interaction matrix ele- 
ments (n\ri2\g 1 713714) by switching to relative and center 
of mass harmonic oscillator quantum numbers, v and N 
respectively. This yields 
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where 4>v{k z ), with </>„ the Fourier transform 

of the z^-th harmonic oscillator eigenfunction. It is eas- 

ily seen that f 2v+l = and f 2v = ( -^f( rj ^) 1/4 \fW- 
The change of basis coefficients are given in Ref. [17] . 
Since (n\n2\Nv) ~ 5N+v,ni+n 2 an d v, v' are even, n\+n2 
must equal, modulo 2, 7i 3 + n 4 to obtain a non-zero inter- 
action matrix element. The 3D contact interaction g can 
be written in terms of the binding energy eb of the two- 
body bound state which always exists in the quasi-2D 
geometry: 
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Here, we simply take N = since Eb is independent of 
the center of mass motion. One can also determine Eb 
as a function of the 3D scattering length a s [12l [13] from 

Eq. (3) using i = ^ (j- - where A is a UV cutoff 
for the 3D momentum that can be sent to infinity at the 
end of the calculation. 

Now if we define the superfluid order parameter 

A qJV=9 ^2 ^"MCq-kna-lCkrut), (4) 
k,ni ,rt2 

and assume fluctuations around this are small, we obtain 
the mean-field Hamiltonian, 
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We further assume that the ground state has a uniform 
order parameter without nodes so that A q jv = <5qo<5jvoAo. 
In this case, Eq. ^ only contains a single unknown pa- 
rameter Aq, so it can be diagonalized to yield 



A 2 

H M F = ^(Ekn-^-^kri) -+ ^ Eknltnvlknv, (6) 

k,n ^ k,n,o" 



where i?kn are the quasiparticle excitation energies. The 
quasiparticle creation and annihilation operators are re- 



spectively given by 

7kr»t = 5Z( Uk "' nC kn't + WnC-kn'4.) (7) 
n' 

7-kr4 = ^( u kri'n c -kn'i — "kn'n4n't) ' (^) 

where the real amplitudes u, u only depend on the mag- 
nitude k = |k| and satisfy J2 n ' (Wkn>n\ 2 + \vk n >n\ 2 ) = 1- 
Note that while they have a well defined spin and mo- 
mentum, they involve a superposition of different har- 
monic oscillator levels. The corresponding BCS wave 
function |*mf) oc Ukna 7kn<r|0), where |0) is the vac- 
uum state for the bare operators Ck na . We then minimize 

A 2 

(H M f) = J2k,n( e kn-fJ>-E kn )--f with respect to A at 
fixed p to obtain the ground state. The value of p is cho- 
sen to keep the density of particles p = 2 ^ k , n |i»k n ' n | 2 , 
and thus the Fermi energy Bp, constant throughout the 
crossover [2"5] . 

For the case where there are only two levels n = 0,1, 
the problem can be solved analytically. Here, Eq. Q is 
greatly simplified since there is no pairing between atoms 
in the n = and n = 1 levels. The quasiparticle disper- 
sions are then E kn = (e kn — p) 2 + (y™Ao) 2 . One can 
now minimize (Hmf) by simply using 8(Hmf) /9Aq = 0. 
Combining this with Eq. |3| yields the implicit equation 
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where we have defined A = A V 00 , the pairing gap m 
the lowest level. Also, the density p = — O(Hmf) /dp, 
and in the regime Ep < u z where Ep = Tzp/m, this gives 
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We see that Eqs. ([£]) and ( 10 ) reduce to the 2D mean-field 
equations |14j in the limit uj z — > oo, as expected. They 
also yield the lowest order correction to the 2D result 
due to confinement (e f /uj z ^ 0): in the BCS regime 
£b/ef <C 1, we have: 
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Thus, A is enhanced by the confinement in this regime 
while p is suppressed. This trend is also observed in the 
full calculation involving many levels (see Figs. [I] and |2]). 

Perturbing away from the 2D limit, we must now in- 
clude multiple levels of the confinement. In general, Hmf 
must be diagonalized numerically to obtain E kn and the 
quasiparticle amplitudes for a given p and A. Equiva- 
lently, one can solve the Bogoliubov-de Gennes equations 
self-consistently, but it is considerably faster to minimize 
the energy (Hmf) directly and it also allows us to take 
into account up to 100 levels. Indeed, we find that higher 
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FIG. 1: (Color online) Chemical potential /j, measured with 
respect to half the binding energy eb for several values of 
the Fermi energy ef/ijJ z . The dashed line is the 2D mean- 
field result |14) . while the dotted curve is the 3D BEC result, 
(i + ~ ^ — 2y ^ F \J Inset: Asymptotic behaviour in 

the BEC regime ef/eb > 1 plotted on a logarithmic scale. 
The solid straight lines are straight line fits to the data with 
gradient —1/3. Error bars for the numerical data are within 
symbol size. 




FIG. 2: (Color online) Behavior of the order parameter 
A throughout the BCS-BEC crossover for different values 
of ef/uj z . The dashed curve is the 2D mean-field result 
A = V2eb£f [H], while the dotted curve is the 3D BEC 

result, A ~ ef yff (^ft 7 ) * ' ^ ne ^ or ^ ne numerical data 
is the same as in Fig. [T] 

harmonic levels are important even for weak interactions 
once Ef/wz shifts away from zero. For the values of eb 
and £f considered in this paper, we can in fact extrapo- 
late the results for fi and A to an infinite number of har- 
monic levels since we find that they both scale linearly 
with the inverse of the number of levels in this limit. 

By incorporating an infinite number of levels, we can 
determine the evolution of the chemical potential /x and 



order parameter A throughout the BCS-BEC crossover in 
a quasi-2D Fermi g cLS, clS depicted in Figs. [T] and [2j The 
lowest density we consider (ef/uj z = 0.1) corresponds 
approximately to the experiments of Ref . [5HS] • Note that 
in the extreme BCS and BEC limits, we must have fj, — > 
Sf and [i — > — £b/2, respectively, a feature which holds 
across all dimensions. In the BCS regime, the behavior is 
in qualitative agreement with the two-level calculation: 
\x is suppressed and A is enhanced with respect to the 2D 
result, with the deviation from 2D being increased with 
increasing ef/^z- However, multiple levels are required 
to correctly capture the dependence on Eb/ef as Ef shifts 
away from zero. We also see in Fig. [I] that the relative 
chemical potential fi+Es/2 exhibits a very steep gradient 
as Eb — > for Ef/u z > 0.5. This illustrates how higher 
levels can lead to a strong deviation from 2D even when 
ef < uj z and eb/ef *C 1 . 

For larger eb/ef, A eventually becomes suppressed 
compared to the 2D result and appears to approach the 
3D mean-field curve in the BEC limit (Fig. [2]). The chem- 
ical potential, however, always remains lower than the 2D 
result, and has a behavior in the BEC regime that is in- 
termediate between 2D and 3D mean field. In the limit 
£b/£f 3> 1 , the pairing gap is no longer given by A and 
the properties of the Bose superfluid are instead encoded 
in /i. In particular, the relative quantity + Eb/% yields 
the mean-field energy for the repulsion between dimers. 
Referring to Fig. [I] (inset), we see that it tends to zero as 
a power law with increasing eb/ef, similarly to 3D and 
in contrast to the 2D mean-field result. This is consistent 
with a dimensional crossover to 3D dimers once Eb ^ u z . 
However, we obtain a power of —1/3 rather than —1/2 
as expected from 3D mean field theory. Indeed, a power 
of —1/2 is also expected for weakly interacting bosons 
confined to quasi-2D since their mean-field energy scales 
as a s /l z ~ ^Juj z /eb [T5]. The discrepancy is likely due 
to the fact that our mean-field approximation does not 
allow for the scattering of dimers in the transverse direc- 
tion since they are constrained to be in the N = center 
of mass mode. The repulsion between dimers has also 
been discussed in the context of a two-channel model for 
the quasi-2D system |19j . 

Deviations from 2D behavior will also be apparent in 
experimental probes of the quasi-2D superfluid. Typi- 
cally, investigations of pairing have exploited RF spec- 
troscopy |20j . where atoms in one hyperfine spin state 
(e.g. I) are transferred via an RF pulse to another hy- 
perfine state that is initially unoccupied. In the ideal sce- 
nario where the final state is non-interacting, the mean- 
field transition rate or RF current is given by 

I RF {u) OC ^ Kn'n| 2 <5(£kn' - M + E kn ~ w) (H) 
k,n' ,n 

where ui is the frequency shift relative to the bare transi- 
tion frequency between hyperfine states. Here, the onset 
frequency of the RF spectrum corresponds to £ , k=o,n=o — 
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ji and is associated with the pairing gap of the superfiuid. 
In the 2D case, ( [IT] ) reduces to I rf {lo) oc ^2<d(uj — Eb) 
and thus the RF pairing gap is simply Eb, as noted by 
Sommer et al. [7J. Perturbing away from 2D, we find 
that the RF spectrum can be substantially modifed by 
the higher confinement levels. In Fig. [3j we see that the 
strongest effects are in the BCS regime, where the pair- 
ing gap is initally enhanced compared to eg. Such an 
enhancement is not surprising given that, in 3D, a pair- 
ing gap can exist even when there is no two-body bound 
state. However, once Ef > w a , the pairing gap drops 
below eb and even becomes negative for small enough 
Eb/ef- This is because the coupling between n = and 
n = 2, and the associated level repulsion (see inset of 
Fig. [3]), reduces the energy E^ n=0 . In this case, the 
lowest energy quasiparticle contains a smaller fraction of 
the n = harmonic level and the RF peak is instead 
dominated by the n — 2 quasiparticle. Thus, the onset 
frequency is no longer an accurate measure of pairing, as 
we can see in Fig. |4j Note that the deep lattices used in 
Ref. [7] correspond to £f/^z ~ 0.03 and thus the pairing 
gap will lie very close to the 2D result, as was observed. 

On the other hand, the experiments in Ref. 8J cor- 
respond to £_f/u; 2 ~ 1.5 and thus the RF spectrum in 
the BCS regime will be strongly affected by the confine- 
ment. In particular, we see in Fig. [4] that the RF peak 
is shifted to higher frequencies compared to the 2D case 
and develops more structure at lower frequencies. Fur- 
thermore, the pairing gap in the BEC regime appears to 
be less sensitive to confinement and closer to the 2D re- 
sult (Fig. [3]) since it is dominated by two-body physics. 
These features are all consistent with the experimental 
observations [8]. We note that finite temperature may 
also play a role in these experiments — see Ref. [21] for 
an alternative explanation based on fcrmionic polarons. 
However, the shift due to confinement appears to be sub- 
stantial at T — 0, with a direction that is consistent with 
experiment, and therefore it cannot be disregarded. In- 
deed, it has also been shown that effects due to confine- 
ment can be significant in spin-polarized quasi-2D Fermi 
gases |2"2"] . 

To conclude, we have constructed a mean-field theory 
for the quasi-2D Fermi gas that is able to capture the 
deviations from 2D behavior resulting from confinement. 
We expect it to provide a benchmark for further inves- 
tigations into quasi-2D Fermi systems. In the future, it 
would be interesting to explore the effects of finite tem- 
perature and the Berezinskii-Kosterlitz-Thouless transi- 
tion in our model. 
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FIG. 3: (Color online) The pairing gap measured using RF 
spectroscopy as a function of eb/ef for a range of differ- 
ent ef/ijJz- In the BCS regime, it deviates substantially 
from the 2D result with increasing ef/^z- In the BEC 
limit, £?k=o,n=o — M must always approach sb, regardless of 
ef/^z- Inset: Lowest quasiparticle dispersions for ef/oj z = 1, 
Eb/ef = 0.1. Note that there is level repulsion between n = 
and n — 2, but no avoided crossing between the n = and 
n — 1 dispersions. 
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FIG. 4: Radio frequency spectra of the quasi-2D Fermi gas 
in the BCS regime with Eb/ef = 0.2 for £f/ui z — 0.5,1.5. 
The RF current Iff is scaled so that the peak value is always 
1. The peak is strongly shifted towards higher frequencies 
compared to the 2D result when ef/^z — 1-5. Note that the 
spectra are calculated assuming that the final state is non- 
interacting. 
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